On the passage from atomistic systems 
to nonlinear elasticity theory 

Julian Brauio and Bernd Schmidllj 
May 31, 2012 

Abstract 

Cn ■ We derive continuum limits of atomistic models in the realm of nonlinear elasticity 

'— ' I theory rigorously as the interatomic distances tend to zero. In particular we obtain an 

ITj ■ integral functional acting on the deformation gradient in the continuum theory which 

depends on the underlying atomistic interaction potentials and the lattice geometry. 
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The interaction potentials to which our theory applies are general finite range mod- 
els on multilattices which in particular can also account for multi-pole interactions 
and bond-angle dependent contributions. Furthermore, we discuss the applicability 
f^ I of the Cauchy-Born rule. Our class of limiting energy densities consists of general 

C^ ■ quasiconvex functions and the class of linearized limiting energies consistent with the 

Cauchy-Born rule consists of general quadratic forms not restricted by the Cauchy 
n , ■ relations. 

^ ■ 1 Introduction 

-4— » 

The main aim of this work is to provide a rigorous derivation of nonlinear elasticity func- 

tionals from atomistic models. The investigation of such discrete-to-continuum limits has 

been an active area of research in continuum mechanics over the last years in particular 

^ I for, but not limited to, elastic interactions. For a recent account on this line of research 
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and a summary of the related literature we refer to the survey article |BLL07) by Blanc, 
LeBris and Lions. 

Classically, the stored energy density in elasticity theory is derived from atomistic 
t~^ ■ models by applying the Cauchy-Born rule: Given a macroscopic deformation y of the 

^^ . elastic body, one assumes that microscopically near every material point x, all the atoms 

deform by just following the macroscopic deformation gradient F = Vy^x). Inserting this 
ansatz into the atomistic potentials then leads to a continuum stored energy density W as 
a function of F € M^^^. Assuming validity of the Cauchy-Born rule, very general and even 
r> ! quantum mechanical interactions have, e.g., been investigated by Blanc, LeBris and Lions 

C^ ■ in |BLL02| . A priori, however, it is not clear if the Cauchy-Born hypothesis does hold true. 

For a two-dimensional mass-spring model, it has been shown by Friesecke and Theil in 
|FT02| that the Cauchy-Born rule does indeed hold true for small strains, while it in general 
fails for large strains. This result has then be generalized to a wider class of discrete models 
and more than two dimensions by Conti, Dolzmann, Kirchheim and Miiller in |CDKM06] . 
A fundamental contribution towards a rigorous derivation of continuum limits in elas- 
ticity has been made by Alicandro and Cicalese in |AC04| . where they prove a general 
integral representation result for continuum limits of atomistic pair interaction potentials. 
It is our main aim, departing from this result, to derive a continuum theory for more gen- 
eral interaction potentials which, in particular, can also incorporate bond-angle dependent 
potentials. Such an extension is desirable in applications, as many atomistic models such 
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as, e.g., the Stillinger- Weber potential, cannot be written as a pure pair potential. In 
fact, the class of potentials our theory applies to is rich enough to model any continuum 
energy density, even if the Cauchy-relations are violated. See below for more details. For 
small strains, such general models have been derived by the second author in |Sch09| . A 
first step in this direction in the nonlinear regime has recently been provided through the 
analysis of a special class of nearest neighbor three point interactions on a two-dimensional 
square lattice by Meunier, Pantz and Raoult, see |MPR11| . and for special interactions 
subordinate to specific triangulations by Alicandro, Cicalese and Gloria [ACGllj . 

Our limiting density will be described in terms of a sequence of cell problems. This is 
related to the homogenization results of Braides |Bra85| and Miiller |Mul87| of nonconvex 
integral functionals. In the setting of a discrete-to-continuum limit for thin films in the 
membrane energy regime, such a limiting cell formula has also been obtained by the second 
author in |Sch08| . Indeed, following the localization method (cf., e.g, |DM93| ). also from 
a technical point of view our reasoning is related to |Bra85| . 

The atomistic systems we will consider are 'generalized mass spring models' as in 
|CDKM06] . For such models, a rigorous simultaneous discrete-to-continuum and nonlinear- 
to- linear limit has been obtained by the second author in |Sch09| . The present paper now 
extends these results in the purely nonlinear setting. While in the linearized regime, as 
shown in |Sch09| . the limiting energy is indeed given by the Cauchy-Born energy, in the 
nonlinear setting this cannot be true in general, cf. |FT02| . Nevertheless, we will obtain 
a definite stored energy density in the limit, which under appropriate conditions will be 
equal to the Cauchy-Born density for small (but finite) strains. 

In order to prove our main representation result we resort to abstract compactness 
properties of F-limits and integral representation results for functionals on Sobolev spaces 
and thus follow the scheme set forth in |AC04| . which is dictated by verifying the hypotheses 
of that abstract approach by the localization method. A few of the arguments in this proof 
can be used with only minor adjustments. There are, however, some major differences 
as compared to the pair interaction case treated by Alicandro and Cicalese. While these 
authors use slicing arguments in order to obtain energy estimates on the usual d x d 
deformation gradients in the direction of interacting pairs, we will have to estimate the 
much higher dimensional d x 2'^ discrete deformation gradients. In fact, as in general our 
discrete energies cannot be recovered by slicing techniques, we will instead work with very 
carefully chosen interpolations of the discrete deformations which encode the full discrete 
gradient on lattice cells. 

More specifically, if C = AZ, is some Bravais lattice, J7 C M a bounded open set with 
Lipschitz boundary that will be viewed as the 'macroscopic' domain occupied by the elastic 
body, whose atoms are at positions eC D 0, we assume that the energy of a deformation 
y : e/^ n O — )■ M'^ can be written in the form 

where x runs over all midpoints of elementary lattice cells of eC inside Q. Here Vy(x) is the 
discrete gradient of y on the corresponding cell Q which encodes all relative displacements 
of atoms lying on the corners of Q. (See Section [2] for precise definitions.) e is the small 
parameter in the system measuring the typical interatomic distance and tending to zero 
eventually in the continuum limit. The rescaling by e is introduced in order to pass from 
units of finite energy per atom to units of finite energy per unit volume. 

In fact, our analysis is not restricted to interactions within unit lattice cells, but also 
applies to general finite range interactions. In such models, the energy is still given as 



the sum over unit lattice cells (£^(ri))°, but the cell energy now depends on the discrete 
deformation gradient Vy(x) G R'^^^ on a larger 'super-cell': 

where the definition of the lattice interior {C^{il.))° is suitably adjusted so that only lattice 
points within may interact. 

Another complication arises when extending our results to general finite-range inter- 
actions on multi-lattices. For such systems, the discrete gradients are augmented with 
additional internal variables describing the relative shifts of the underlying single lattice. 
With the help of a mixed Sobolev/Lebesgue space representation theorem we are then 
led to a boundary value/mean value cell formula for the limiting energy density. This 
cell formula is in fact related to the cell formula derived in |Sch08| for thin membranes 
where internal variables measure relative shifts of the thin film's layers. On multi-lattices 
e({0, si, . . . , Sm} + >C) with m shift vectors si, . . . ,Sm € K'' the general discrete energy 
functional then reads 

Fs{y,S,n)=£'^ ^ T^super-cell(Vy(2;),s(x)), 

xe(C',{n))° 

where y : sCCiQ, ^ W^, s : eCCiQ, ^ W^^"^. For notational convenience we will restrict to 
simple unit cell interactions on a Bravais lattice for the largest part of the paper and only 
comment on the necessary modifications in the more general case at the end of Section \E\ 
Our main results are summarized in the following theorems. The necessary assumptions 
Assumptions [H [2] and [3] on the cell energy are specified in Section [2j (Assumptions [1] and 
[2] are nothing but standard p-growth assumptions on Wcdi-) 

Theorem 1.1 (F-convergence) . Suppose Assumptions {I\ and\^ are true. Then Fi,{-.,Q) 
F(LP(ri;M"'))- andT{Ul^^{9.;W^)/M.)-converges to the functional F , defined by 



F{y) = { 



W^cont(Vy(x))dx, ify(^W^-P{n-] 
n 
, oo otherwise, 

where the continuum density VFcont : M — > [0, oo) is given in terms ofWcdi by 
W,ont{M) = -^ ^lirn^ i^ ^^M Yl VFceii(V2/(x)) : y € Bi{A{0, N)^, yu) 

Here Bi{A{^, N) ,yM) is the space of lattice deformations of Ci H A[0, N) with linear 
boundary conditions M on dCi{A{0, 1)'^), cf. Section 2] 

As in non-convex homogenization (see |Miil87| ). in the representation result for VFcont 
it is nesessary to minimize VFceii over larger and larger cubes and the limit is in gen- 
eral not obtained for finite N. A simple 2d example for this effect is given by a square 
lattice where nearest neighbor atoms interact via a harmonic spring potential: Fi,{y) = 
2 S|x-x'|=e(ly(''^) ~ y{^')\ ~ ^)^ (which can be written in the above form). This is a simpli- 
fied version of the example discussed in |Sch08| Sect. 4.4]. The arguments sketched there, 
which amount to considering deformations 

yixi,X2) = M{xi,X2) + cri{xi) + 0-2(x2) 



for compressive boundary conditions, with suitable 2-periodic functions cji and (T2 like 






i = 1, 2, and suitably modified on the boundary, show that 

^infj J2 W,,n{Vy{x)):yeBi{{0,N)'',yM] 

[xG(£;((0,Af)d))° 

converges to (max{0, \m.i\ — 1})^ + (max{0, \m.2\ — 1})^ with error bound 0{N^^), where 
m.j denotes the jth column of M. Evaluating, however, the energy of the bonds on and 
close to the boundary it is easily seen that, for compressive boundary conditions |m.i[ < 1 
or \m.2\ < 1, this limiting energy is always over-estimated by a constant times N~^. 

Theorem 1.2 (Compactness). Under the assumptions of Theorem ] 1.1\ if y^ is a sequence 
with equibounded energies Fi;{yfr, Vi) and Vt is connected, then there exist a sequence e^ — )■ 
and y G W^^P{n;R'^) such that y^^ -^ y in L^^^{n;R'^)/R. 

Of course, if Q is not connected, one has compactness in L^^^ up to translation on every 
connected component. 

Analogous results hold true under boundary conditions g G M^^'°°(R'^;R'^). Let F^ and 
Fi denote the functionals obtained from F and Fg, respectively, with values set to infinity 
if the boundary conditions are not met, cf. Section |4] 

Theorem 1.3 (F-convergence) . If Assumptions\^and\^are true, thenFi{-,fl) r(L^(J7;M'*))- 
converges to the functional F^ . 

Theorem 1.4 (Compactness). Under the assumptions of Theorem \1.3\ if y^ is a se- 
quence with equibounded energies Fi{yg,Q,), then there exist a sequence Sk ^- and 
y € W'^^P{n;R'^) with y = g on 89. such that y^^ -^ y m U'{9.;R'^). 

A standard argument then yields that almost minimizers of Fi{-,Vt) converge to min- 
imizers of F^{-,Q) and almost minimizers of F^{-,Q) up to translation converge to mini- 
mizers of F{-,n), more precisely: 

Corollary 1.5 (Convergence of almost minimizers). Suppose Assumptions U\ and\M are 
true. Then every sequence of almost minimizers of Fi;{-,0,) for connected il. is compact 
in L^^^{Q;R'^)/R and every limit is a minimizer of F, while every sequence of almost 
minimizers of Fi[-,Q,) is compact in L'''{Q,;R ) and every limit is a minimizer of F^ . 

It is not hard to include body forces in the energy expression as these will only be 
continuous perturbations of the energy functional which converge uniformly on bounded 
sets and thus preserve F-convergence by general theory. 

We also remark that the point why the theory can be adapted to the case of general 
finite range interactions, is that in this case Wcclh while naturally still bounded from above 
by the discrete gradient through Assumption [J] from below has to be bounded only in 
terms of the discrete gradient on the unit cell. See Section [S] for details. For general finite 
range interactions on multi-lattices we state here only the analogue of Theroem II. H as 
in the cell formula there are now additional internal variables that need to be taken into 
account. Theorems 11.21 [l~3] and \TM and Corollary 11.51 extend in a straightforward manner. 



Theorem 1.6 (F-convergence) . Suppose VKjuper-cell satisfies the growth assumptions (as 
stated in Section\^. Then Fe{-,-,n) r{LP{n;R'^) x w-Li{n;R'^''"')) -converges to the func- 
tional F , defined by 



F[y,s) = { 



n 
,00 otherwise, 

where the continuum density Wcont : K'^^^ x R^,^"^ -^ [0, oo) is given in terms o/ Wsupcr-cell 
by 

Wcont{M, So) = rj^^ j^^J^l^ iV^ ^"M ^ VFsuper-cell(Vy(x), s(x)) : 

' ^ ' ^°° I xe{C[{A{0,N)'i))° 

{y,s)eBi{A{0,Ny,yM,so)\. 

Here Bi{A{0, N) , yM , sq) is the space of lattice deformations of Ci fl A{0,N) with 
hnear boundary conditions M on dCi{A{0, 1)'^) for y and average sq for s, cf. Section [5j 

As macroscopic deformations are solely given in terms of a deformation mapping y € 
W^'P{Q,; W^), we are also interested in the effective macroscopic energy density obtained by 
minimizing out the internal variables s: 

Theorem 1.7. For every y G W^'P{Q;W^) we have 

min Wcont(yy{x),s{x))dx = / min Wcont{'^y(.x),s) dx. 

n n 

Moreover, F/~™'"(-, J]) := inf Fe{-,s,n) T{LP{n;M.'^))-converges to the functional F""'^''' , 
defined by 



F"-™°(y) 



[ min Wcont{'^y{x),s)dx, if y £ W^'P{n;R''), 
n 
, oo otherwise. 



Returning to our basic setting on a Bravais lattice, under an additional assumption, we 
can calculate the limiting density for small strains explicitly by the Cauchy-Born rule: 

Theorem 1.8. In addition to Assumptions[l\andl^ suppose that WccW satisfies Assumption 
O Then there is a neighborhood U of S 0(d) such that Wcont is given by 

W,ont{M) = Wcb{M) := —L^WcelliMZ). 

\aetA\ 

for all M eU. 

Here Z G R"'^^ is a 'discrete identity matrix', see Section [2] for details. 

As Wcont arises as the energy density of a F-limit it has to be quasiconvex (cf . Section [2] 
for the definition of these concepts) . The next proposition shows that our class of atomistic 
interactions is rich enough to model any quasiconvex energy density in the continuum limit. 



Proposition 1.9. Suppose V : M.'^^'^ — )■ M is quasiconvex with standard p-growth 

c \M\P -d < V{M) < c"{\Mf + 1) 

for some constants c,c',c" > and all M € R . Then there exists a cell energy Wceii 
satisfying As sumptions {1\ andl^ such that Wcont = V- 

We remark that, by way of contrast, a restriction to pair interaction models will only 
lead to a restricted class of limiting continuum energies, as can be quantified in terms of 
the so-called Cauchy relations: If the Cauchy-Born rule applies (e.g., due to Assumption 
[3|), an atomistic interaction energy 









yields the continuum density 



^cb(m)=--i^5:ii.. 



Idet^l 






Assuming Vi^i is smooth and, for large \x\, sufficiently rapidly decreasing a direct calculation 

yields 

d 

D'^WcB{ld){M,M)= Y^ Cijkimijruki, 

i,j,k,l=l 

where the elastic constants CijM are given by 

^ \~~^ -r r// n i\XiXjXjsXl . I XjXlOki XiXjX^Xl \ 

While the symmetry relations Cijki = CkUj and Cijki = Cjiti naturally follow from the 
symmetry of the Hessian Z)^VFcB(Id) and frame indifference of VVcB) the particular form 
of WcB in addition gives djki = cukj = Ckju for every i,j, k, I. 

In the 3-dimensional setting of elasticity theory these additional relations lower the 
dimension of admissible elasticity tensors from 21 to 15 (symmetric in all indices) and so 
can be written as 6 equations, the Cauchy-relations, namely 

C1122 = C1212, C2233 = C2323 5 C3311 = C3131, 
C1I23 = C1213) C223I = C232I, C3312 = C3132. 

The question whether in elasticity theory the Cauchy-relations hold true (rari-constant 
theory) or fail (multi-constant theory) had been under discussion for quite some time in 
physics and was finally decided by experimental data in favour of the multi-constant theory 
(for some experimental data and further physical considerations cf. |Hau67| ). This means, 
in particular, that the interaction in a lattice is a complex multibody interaction which 
cannot be reduced to pair-potentials. Our model in this paper using general cell energies 
is not limited by the Cauchy-relations: 

Proposition 1.10. Suppose Q : M ^ R is a positive semidefinite quadratic form which 
is positive definite on the symmetric dxd matrices and vanishes on antisymmetric matrices. 
Then there exists Wceii satisfying AssumptionsUl andl^ such that 

lD^WcB{ld){M,M) = Q{M). 



The paper is organized as follows. In Section [2] we first introduce the discrete model 
and review some basic facts on T-convergence and integral representations of functionals 
on Sobolev spaces. In Section [3] we then proceed to state precisely and prove a general F- 
compactness and representation theorem. This in particular requires a number of technical 
preliminaries in order to investigate discrete deformations. A version of this representation 
result for boundary value problems is then provided in Section HI Finally, the limiting 
stored energy function is identified in Section [5] through minimizing a sequence of cell 
problems, leading to the main discrete-to-continuum convergence result and the proofs of 
the results stated above. 

2 The model and general preliminaries 

In this section we introduce the atomistic model and recall some general facts on F- 
convergence and integral representation results required by the localization method. 

2.1 The atomistic model 

Let C CM. he a Bravais lattice, i.e., there are linearly independent vectors vi, . . . ,V(i such 
that 

C = {nivi H h UdVd \ni,. . . ,nd &Z} = AZ"^, 

if we set A = {vi, . . . ,Vii)- The scaled lattices C^ = sL partition M into the e-cells 
z -f ^[0, e)"^ (2: G £e). Let Qe{x) denote the e-cell containing x. The centers of the cells 
are C^ = C^ + A{-^, ■ ■ ■ j 2) ^^^ ^^ denote by x the center of the cell containing x. These 
centers give a convenient labeling of the cells. Furthermore let zi, . . . , Z2d be the points in 

A {-i, i}^ and Z:={zi,...,Z2d)e R'^''^' . 

For a set U CM.'^ we define the following lattice subsets in the spirit of its closed hull, 
interior or boundary with respect to eC or its corners eC by 



4(C/) = {x e 4 I Qe{x) n [/ / 0}, CeiU) = C',{U) + e{zi,. . . , z^,], 



{C',{U)r = {x e 4 I Qe{x) C f/}, {Ce{U)y = {C',{U)r + e{zi, . . . , Z24, 

94(c/) = 4(t/)\(4(c/))°, dc,{u) = dc',{u) + 6{zi, . . . , z^,]. 

Furthermore let 

U'= [j Qs{x), Ue= U ^^(^)- 

xec',(u) xe(a(u))° 

A lattice deformation should be thought of as a mapping /^^ fl $7 — t- R'^. Choosing a 
suitable extension (e.g., by 0) and piecewise constant interpolation, we can and will assume 
that the lattice deformations Se(r2) are the functions O — )■ W^, which are constant on every 
cell Qs{x), X G £^(J7). (This will not change the energy, see below.) 

If we have a deformation y € Bs{0,) and x € H.^, we set yi{x) = y{x + £Zi), 

r\d 
1 - 1 

yi^) = ^^yii^) and Vy{x) = -{yi{x) -y{x),...,y2d{x) -y{x)). 
1=1 

Let A{U) be the set of all bounded open subsets oi U C W^ and Al{U) the set of all 
those, that have a Lipschitz boundary. In the following, we will consider a set Q G Al{R ) 



and the energies F^ : LP{i^; W^) x A{0,) — ?• [0, oo] for some fixed 1 < p < oo, defined by 



Feiy,U) 



oo otherwise. 



(2.1) 



In this definition the energy only depends on the values of y in (£.^{U))° C C^ CiU. Of 
course, there can be some points in C/^nU which we do not use at all, but this is negligible 
if we impose Dirichlet boundary conditions as we will do later on. 

We make some assumptions on the cell energy Wcdi: M 
discrete gradient can take values precisely in the space 



[0,oo). Note, that a 



Vn 



F G M'^'^^ : Y^ aij = 0, for every i = 1, 



Therefore, we are only interested in the values of VFceii on Vq. 
Assumption 1. There are c, c' > such that for every F € Vq 

Wccii{F)>c\Ff-c'. 
Assumption 2. There is a c > such that for every F € Vq 

^cell(i^) < c{\Ff + 1). 

While these conditions are supposed to hold true for all our results, we also state a 
third assumption, which, if satisfied, allows for an application of the Cauchy-Born rule 
locally near SO{d). The so-called Cauchy-Born energy density is defined by letting each 
atom follow the macroscopic gradient: 



Wcb{M) :-- 



1 



Idet^ 



■W,,iiiMZ) 



for M G 



pdxd 



Assumption 3. (i) Wccii : 



0^x2" 



is invariant under translations and rotations, i.e. 



ix2'' 



for F € 

W,,n{RF + ic,...,c)) = W,en{F) 
for ah R G SO{d), c G M*^. 
(ii) Wccii{F) is minimal (= 0) if and only if there exists R G SO{d) and c G M such that 

F = RZ + {c,...,c). 

(iii) VFccii is C^ in a neighborhood of SO{d) := SO{d)Z and the Hessian D'^Wcc\\{Z) at the 
identity is positive definite on the orthogonal complement of the subspace spanned 
by translations (c, . . . , c) and infinitesimal rotations FZ, with F-^ = —F. 

(iv) p > d, which together with Assumption [1] implies in particular that Wceii satisfies the 
growth assumption 

Fi->oo \Fr 
Fevo 



2.2 F-convergence and integral representations 

In our analysis, we consider energies on discrete systems depending on a small parameter e, 
the scale of the lattice spacing. To make the limit for e — )• precise and gain some knowl- 
edge about the behavior of associated minimizers, we will use De Giorgi's T-convergence. 
We recall the definition and some basic properties that will be needed in the sequel. 

Definition 2.1. Let X be a metric space and F„, F : X ^ M = MU { — oo, oo}. We say F„ 
r(X)-converges to F {Fn ^ F), if 

(i) (liminf-inequality) For every y,yn €z X with y„ — t- y, we have 

F(y)<liminfF„(y„), 

(ii) (recovery sequence) For every y G X, there is a sequence yn € X such that 

F{y) > limsupF„(y„). 

If (-F£)e>o is a family of functionals depending on a positive real parameter e, we say 
Fg r(X)-converges to F, if for every sequence e„, > converging to 0, we have F^^ — > F. 
We will also use the F-limsup and the F-liminf, given by 

F'{y) = r(X)-liminf F„(y) = inf{liminf F„(y„) : 2/„ ^ y in X}, 

n— >oo n— >oo 

F"{y) = r(X)-limsupF„(y) = inf{limsupF„(y„) : ?/„ -> y in X}. 

71— ^-oo n— ^oo 

Note that (i) is equivalent to F < F' and (ii) is equivalent to F > F" . Hence, Fn — > F if 
and only if F' = F" = F. Furthermore, we see that F-convergence is a pointwise property, 
so we can speak about F-convergence at a specific point. 

In the following proposition we assemble some basic properties of F-convergence that 
we will not prove here. 

Proposition 2.2. (i) The infima in the definitions of F' and F" are actually attained 
minima in M; 

(ii) every sequence of functionals on a separable metric space, like LP{U;M. ), has a F- 
convergent subsequence; 

(Hi) F' , F" and F are lower semicontinuous with respect to convergence in X . 

(iv) T -convergence satisfies the Urysohn property, i.e., Fn V -converges to F, if and only 
if every subsequence of Fn has a further subsequence, that T-converges to F; 

(v) if Fn T-converges to F and Gn converges uniformly on bounded sets to a continuous 
functional G, then Fn + G„ T-converges to F -\- G. 

In view of applications, the most interesting property of F-convergence is the following 
theorem. 

Theorem 2.3. // Fn T-converges to F and sequences yn in X with equibounded Fn{yn) 
are pre- compact then F attains its minimum on X and we have 

umiF{x) = lim inf Fn{x). 



Furthermore, let y„, (^ X be a sequence with 

FniVn) -^ lim inf F„,(x), 

n— >ooxeX 

then the limit of every converging subsequence of y„ is a minimizer of F. 

For proofs of Proposition 12.21 and Theorem 12.31 see, e.g., |DM93| . 

Returning to our specific setting, for a sequence e„ > such that e„ -^ 0, we define 



F'{y, U) := T{LP{^- W''))- Urn inf F,„ (y, U) 

n— >oo 

= min{liminf Fe„(y„, U):yn^y in L^i^; M'^)}, 

n— >oo 

F"{y, U) := r(LP(0; M'^))- lim sup F,„ (y, C/) . 

= min{nmsupFs„(y„,C/): y„ -> y in LP(J7;M'^)}. 

n— )-oo 

The Hmiting functionals we will encounter in the next section are integral functionals 
of the form 

/: W^i'P([/;]R'=) ^ [0,oo], I{y) = [ f{Vy{x))dx 



u 

with 1 < p < oo,U G A{W^), f : M'^^'' — )■ [0, cxd) continuous. Recall that a Borel measurable 
and locally bounded function / : M — >■ R is quasiconvex, if 



f{M)<jf{M + V^{x))dx, 



u 

for every nonempty U € ^(M''), M € M''^'^ and (/? G Wq''^{U;R''). In our analysis, the 
quasiconvexity of / will be due to the following result. 

Theorem 2.4. /// is sequentially weakly lower semicontinuous in VF"'^'^(C/;M ), then f is 
quasiconvex. 

A detailed discussion of quasiconvexity and related properties, including proofs of the 
above statements, is given, e.g., in |Dac08| . 

In order to guarantee that indeed our limiting functional is an integral functional, we 
will resort to the following general integral representation result on Sobolev spaces. 

Theorem 2.5. Let 1 < p < oo and let F: W^^P{n;R'^) x A{Q) -^ [0, oo] satisfy the 
following conditions: 

(i) (locality) F{y, U) = F{v, U), if y{x) = v{x) for a.e. x € U; 

(a) (measure property) F{y, •) is the restriction of a Borel measure to A{'^); 

(Hi) (growth condition) there exists c > such that 



F{y,U) <c j'\Vy{x)\'P + ldx; 
u 

(iv) (translation invariance in y) F{y, U) = F{y + a, U) for every a S M'^ ; 
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(v) (lower semicontinuity) F{-,U) is sequentially lower semicontinuous with respect to 
weak convergence in W^''^{Q;M. ); 

(vi) (translation invariance in x) With yuix) = Mx we have 

F{yM,Br{x)) = F{yM,Br{x')) 

for every M G M"'^'', x,x' e^ and r > such that Br{x),Br{x') C O. 

Then there exists a continuous f: M — )• [0, oo) such that 

< /(M) < C(l + \Mf) for every M G R"^""^ and 

F{y,U)= f f{yy{x))dx. 



u 

A proof can be found in |BD98| pp. 77-81] or in the scalar-valued setting, which is 
essentially the same, in |DM931 pp.215-220]. 

To show the measure property in the previous theorem, we will use the following lemma. 

Lemma 2.6 (De Giorgi-Letta). Let X he a metric space with open sets r. Assume that 
p: T — )• [0, oo] is an increasing set function such that 

(i) p(0) = 0, 

(ii) (suhadditivity) p{U UV) < p{U) + p{V) for all U,V er, 

(Hi) (inner regularity) p{U) = sup{p{V) : V (z t,V CC U} for all U (z t, 

(iv) (superadditivity) p{U UV)> p{U) + p{V) for all U,V e t with UnV = $. 

Then the extension p of p to all subsets of X , defined by 

p{E) = inf{/9([/) :UeT,EcU}, 

is an outer measure and every Borel set is p-measurahle. 

For a proof see, e.g., |FL07l pp. 32-34]. 

3 A general representation result 

In this section we will prove a general compactness and representation result for sequences 
of discrete deformations. For pair interactions, the following theorem has first been estab- 
lished by Alicandro and Cicalese in |AC04| . 

3.1 Statement of the representation result 

Theorem 3.1 (compactness and integral representation). Suppose Assumptions\^ and\^ 
are true. For every sequence e„ > such that e„ — )> 0, there exists a subsequence £„,, 
and a functional F: LP{Q;W^) x A{0,) — )■ [0, oo] such that for every U G A{0,) and y G 
W^'P{n;R'^) the functional Fe„^{-,U) r{LP{n;R'^))-converge to F{-,U) aty. Furthermore 
there exists a quasiconvex function f: R'^^'^ — > [0, oo) satisfying 

c \Mf -c' < f{M) < c'(|M|P + 1) 
11 



for some c,c' > such that 

F{y,U) = J f{Vy{x))dx if y e W^'P{n;R''). 
u 

In addition, ifUG ^/,(0) (in particular, ifU = 0,), we have 



F{y,U) 



, ^^v j, 



U 

l,CXD 



otherwise, 
and the junctionals Fg,^ (•, U) T-converge to F{-, U). 

3.2 Interpolation 

We now define the continuous and piecewise afiine interpolation y of y, similar to |Scli09| : 

First consider the cell A \—2i 2] ^'^'^ V- ^ {~2' 2} ~^ ^'^^ ^"^ every 0-dimensional face of 
the cell just take y = y. Now assume we already have chosen a simplicial decomposition on 
every [k — l)-dimensional face and have interpolated affine there. Let F = cojzj^ , . . . , Zj j. } 
be a A;-dimensional face. Set 

2^ 



^X]^^-' y(^) = ^Z]2/(^i™)- 



m=l 



m=l 



To complete the induction, we decompose F into the simplices co{«;i, . . . ,1/;^,^}, where 
co{wi, . . . ,Wk}\s a simplex belonging to a simplicial decomposition of an (n— l)-dimensional 
face. Define y to be the interpolation affine on every constructed simplex. If y € j8e(r2), 
we get y on O^ by interpolating as above on every cell. 

The following proposition is about the relation of Vy and Vy. 



Proposition 3.2. There are C, c > such that for every x G O^ and y € Be{^) 
c|Vy(a;)|^< i \Vy{x')f dx < C\Vy{x)f . 



(3.1) 



Proof. Every Zi belongs to some simplex K of the construction. Choose a G K°, where 
the gradient is well-defined. Since the interpolation is linear on K, we see that 



y(x + ezi) - y 



= \\/y{a)z,f 
<C\Vyia)r 
= Cj^J\Vy{x')fdx' 

K 

<C^QlM I IVyix')]" dx' 
\K\ J I ^^ ^1 

Qe{x) 

<C i \Vy{x')\^ dx' , 

Qeix) 
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where C is independent of x, e and y. We immediately get the first inequahty. For the 
second inequahty we prove by induction over k that for every fc-dimensional simplex S = 
co{z, Zi,wi, . . . , Wk~i} in the construction regarding Qe{x) we have 

|Vy(a)Pyr<C|Vy(x)|^ (3.2) 

for every a & S, where Py is the projection on V = span {z — Zi,wi — Zi, . . . , Wk~i — Zi}. 
The case A; = 1 is clear since then for some j we have V = span {zj — Zi} and 

Vy{a){zj - Zi) = Vy{x){ej - ej). 

If the statement is true for k — 1, we immediately have (|3.2|) for 

V' = span {wi - Zi,. . . , Wk-i - Zi}. 

But as in the k = 1 case we also have (I3.2p for V" = span{z — Zi} = span {zj — Zi}. Let 
us define 

1 1 7) 1 1 U/l -I- U/'l 

II lly — I I ^ I I ' 
if V & V , v' € V' and v" € V" such that v = v' + v" . This is a norm on V and hence we 
can calculate using the equivalence of all norms on finite dimensional spaces 

\Vy{a)Pvf < Csup{\Vy{a)vf : v e V, \\v\\y < 1} 

< C(sup{|Vy(a)i;'|^ : v' G V , \v'\ < 1} + sup{\V y{a)v'\^ : v G V\ \v'\ < 1}) 
<C|Vy(x)|^. 

Since we have only finite many possibilities for V, V, V", this C can be chosen independent 
of them, which concludes the induction. Take k = d and integrate to get the result. D 

Proposition 3.3. Let En > 0, with e„ — )■ 0, y,i G Be^{^) and y G LP{Q.]W^) such that 
Vn ^ y in LP(r2;R'^). For every V CC il, we then have yn ^- y in LP{y;W^). 

Proof. It is enough to show ||y„ ~ ynWiviv-R'^) ~^ '-'■ ^^^ Aj : M — )■ [0,1] denote the cell- 
periodic functions such that 

yn{x) =^^^i\ — ] yn{x + EnZi) 

= ^ Aj ( — I yn{x + en{Zi - Zi)), 

i=l ^^"^ 

where without loss of generality we have chosen a numbering of A{— 2, 2} such that 
Z\ = ^(—2, . . . , — 2)- Of course, A^ > and the Aj add up to 1 in any point and so for n 
large enough 

V V V-^ " / 

< / max \yn{x) - yn{x + Enizi - zi))\^ dx 

J i=l,...,2'* 
V 

- X] / ly'«(^) ~ yn{x^En{zi - Zi))f dx. 



-w 
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But the last term goes to since for every i € {1, . . . , 2"^} 

WVn - yn{- + en{Zi - ^l)) |lLP(y;Rd) < 2 ||y„ - y||ip(Q.ud) + Wv - Vi' + Sn{zi -^l))||ip(y. 

^0. 



D 



3.3 Preliminary lemmata 

We proceed to collect further lemmata. We will use them later to prove the requirements 
of Theorem 12.51 In the following, fix some sequence of positive real numbers e„ — )> 0. 

Lemma 3.4. Suppose AssumptionU\is true. If y € U\Q.;'M. ) and U S A{^) are such that 
F'{y,U) < oo, then y G W'^'P{U;R'^) and 

F'{y,U)>c\\Vy\\l^^^,^^,.,,^-c'\U\, (3.3) 

for some c, c' > independent of y and U . 



Proof. Let y^ ^^ y vn L^(il;]R ) such that liminf F£^(y„, ?7) < oo. For some subsequence 
7ifc, we have 

lim Pen, {ynk,U) = liminf Fe„(y„, U), 

yrik S Bs„ (U) and F^^ {yn,^,U) < M < oo for some fixed M > 0. By Proposition HTHl 
we have y„j, — )> y in LP(V,M.) for every V CC U. Furthermore, by Assumption [T] and 
Proposition 13. 2| we get 



M > Fe„^ {yn,,U) = ei^ ^ TVceii(Vy„,(x)) 



^^i^L,(u)r 



>4 E (HVyn,(x)r-c') 



^^i^L,(u)r 






^6(^U,(C/))° 




j \Vyn,{x')\''dx'-c'\ 



We thus obtain 



/ \\/yn,ix')f dx' < M + c' \U\ 



c 

hence the gradients are bounded in LP(V;W^). By the properties of weak convergence 
on Sobolev spaces this means y € W^'P{V,R'^) and Vy„j. ^ Vy in LP{V;M.'^). Weak 
sequentially lower semicontinuity of the norm yields 

c l|Vy||^j,(^.igdxd) < i™™f ^£n(yn, U) + c \u\ , 
but the right hand side is independent of V, thus y € VF^'^(C/, M'^) and 

c II Vy||^,(^^^,,,^ < li^minf F,„(y„, [/) + c' \U\ . 
The definition of the P-lim inf now yields the lemma. D 
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Lemma 3.5. Suppose Assumption\^ is true. Then there is a C > such that for every 
V G Al{^), U G A{V) and y G LP{n;M.'^) n W^'P{V;M.'^) we have 



F"{y,U)<C{\\Vy\\l^^^,^^,^,^ + \U\ 



(3.4) 
Proof. We first prove ([0|) for every y G C^(]R'';M'^). For x G C'^^ and a G Qeni^) define 

yn{a) = y{x). 

Thus y„ G Be„{U) and since y is uniformly continuous, we have y-n ^ U uniformly and 
hence in L^(f];]R ). By Taylor expansion we have 



yn{x) - yn{x + EnZi) 



y{x) - y{x + EnZi) 



<C{\Vy{x)\+en\\V^y\ 



With Assumption [2] we can calculate 



F,Syn.U) = ei Y. ^ceii(Vy„(x)) 

<Cei Y. (|Vyn(x)|^ + l) 

< C \U\ + Cei Y (iVy(^)l' + 'n II V^yf^) 

< C \U\ + C" \U\ el \\vh\L + Cei E 1^2/(^)1' 

Furthermore for every x' G Qe„{x), x G (£^^(C/))° 

|Vy(x)r < C(|Vy(x')|^ + \Vy{x) - Vy{x')\'') 
<C{\Vy{x')\'' + el\\V^y\\''j, 

and, by integrating over x' and summing over x, we get 



'71 / J 

x€iaju)y 



\Vyix)\' <cil \Vy{x')\' dx' + \U,^\eI ||V\||^ 

\^.„ / 

<C ij \Vy{x')\'' dx' + \U\el\\V^\A . 



Putting the two inequalities together and letting n ^^ oo, we obtain 



limsupF,„(y„,C/)<C ||Vy||^ 



LP(C/;R''x'') 



+ \u\ 



So by the definition of the F-lim sup we have ()3.4p . 

The general case follows easily: Since V has Lipschitz boundary, we can take yk G 

C^(]R'^;M'^) such that yk ^ y in W^''p{V;M.'^). Then we have by lower semicontinuity of 
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F%.,U) 



F"{y,U) <limmi F"{yk,U) 

< i™™f c* (llvyfcllip{c/;Rdxd) + 1^1 

D 

Lemma 3.6. Suppose Assumptions {I\ andl^ are true. Let U,V,U' € A{il.) be such that 
U' CC U. Then for every y e VFi'P(f^;M'^) 

F"{y, U'UV)< F"{y, U) + F"{y, V). 

Proof. Without loss of generality, we can assume the terms on the right hand side to be 
finite. According to the properties of the F-limsup it is possible to find sequences Un,Vn 
such that 

limsupF,„(n„,C/) = F"(y,C/) 

n—^oo 

limsupFe„{vn,V) = F"{y,V) 

n— >-oo 

Vn^yin LP{n;R'^) 

For n large enough F^^^Un, U) and F^^^Vn, V) are bounded and Un S Be„{U), Vn S B£^{V). 
Fix TV G N, iV > 5 and then define D = dist(C/', U") and Uj = {x e U : dist(x, U') < 
^}. Choose cut-off functions ipj such that 

(Pj{x) = 1 \/x £ Uj, 

¥., ecr(f/i+i;[0,l]), 

„^ n 2iV 

V(i9,- < . 

II r J lloo — jj 

Next we define 

Wn,j{x) = ipj{x)Unix) + (1 - ipj{x))Vn{x) 

and calculate 

Wnj{x + EnZi) - Wnj{x) _ _ ^_____^U„(3; + EnZi) - Un{x) 



ipj{x + enZi 



£ 



n 



, /I / ; \^Vn[X -\- £nZi) Vn\X) /n rr\ 

+ {l-ipj{x + £nZi)) (3.5) 

+ [Unix) - Vnix))^ '- . 

£n 

To estimate F^^{wn,j,U' U V), we have to look at (C^ (U' U V))° . Clearly, if x is in 
{C'^^{Uj))°, then Vwn,j{x) = Vun{x) and if x is in {/:'^^{V\Uj+i))° , then Vwnjix) = 
\/Vn{x). To control the other cases, observe that for n large enough diam((5£„) < jj^ and 
thus 

(4„(c/'uy)r c (4„(c/,)r u(4„(nt^)r u(4jyn(c/,+2\i^))r 
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for every j € {2, . . . , A^ — 3} and n large enough. With Wj = V (1 {Uj^2\Uj-i), we then 
have 

xeiC',JU'uv)r 

< Fe„{Un, U) + F,„{Vn, V) + 8^ ^ W,,n{^Wnj{x)) . 

' V ' 

We now have to estimate Sj^n- For all n large enough, use first Assumption [2] and then 
(|33]) to get 

Sj,n<Cei Yl (|Vzz;n,i(2;)r + l) 

xe(C',Jw,)r 
< Cei Y (|Vn"(^)r + l^Mx)]" + K{x) - Vn{x)rNP + 1) 

xGic',jw,)r 



<C \VUn{x)f+\VVn{x)f + NP\Unix)-Vn{x)f + ldx, 

because of the gradient of ip being bounded by CN and Proposition 13.21 Averaging over j, 



we get 

Ar-3 
i 

N 



Ar-3 

3^ Y Sj,n < Cj;^^ I |V^n(x)r + |Vi)„(x)|P + ldx + NP I \Un{x) - 7;„(x)|P dx. 

(3.6) 
Of course we can always find a number jij^) such that 

.. Af-3 

i=2 

By Proposition 13.21 and Assumption [H the first integral in (13. 6p is bounded, but 

for n — > oo, hence 

limsup5j(.„)^„ < 



N -\ 



If we define y„ = Wn.jin)-, then obviously y„, G Bs„{U' U y) and y-n ^ V in LP(ri;M ). We 
have 

F"(y, [/' U y) < limsupF,„(y„, U' U F) 

n— >oo 

< hm sup Fe„{un,U) +\imsupFe„{vn,V) +limsup5j(„)^„ 

n—s-oo n ' 

<F"{y,U) + F"{y,V) + 



n—^oo 

c 



N -A 
Letting A^ — ?> oo, we get the conclusion. D 
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Lemma 3.7. Suppose Assumptions {J\ and\^ are true. Then for every V € Al{^), U € 
A{V) and y e LP(0; W^) n W^^p{V; W^) 

F"{y,U)= sup F"{y,U'). 
u'ccu 

Proof. Since F"{y, •) is an increasing set function, we only have to show '<'. 
Let (5 > 0. Then take a U'" CC U such that 

Choosing U', U" such that 

U'" CC U" CC U' CC U, 

we can calculate 

F"{y,U)<F"{y,U"^U\lj^) 

< F"{y, U') + F"{y, U\W) 

< F"{y, U') + 6C, 

where we used Lemma 13.61 and Lemma 13.51 D 

Lemma 3.8. Suppose Assumptions {J\ and\^ are true. Then for every V € Al{^), U G 
A{V) and u,v e LP{Q;R''-) n VF^'P(y;M'') such that u{x) = v{x) for almost every x eU, 
we have 

F"{u,U) = F"{v,U). 

Proof li u = V a.e. in U then for U' CC U we have F"{u, U') = F"{v, U'). To see this, 
just change any approximating discrete sequence of u outside of {U'Y'^ such that the new 
sequence converges to v. But this is enough by Lemma 13.71 D 

3.4 Proof of the representation result 

Now, we can finally prove the compactness result: 

Proof of Theorem VJ.l\ First we find by a suitable diagonal argument a subsequence Fg,^ 
such that we get F-convergence for every U G A{VL). For this we define 

Ai = \u d^: C/=Qfi,,(xi),x, GQ^r, GQ,r, >0,iV€N. I 

The set Ai is countable and we can write ^i = {C/i, C/2) • • • }■ Now choose subsequences 
as follows: 

Ff,^{-^Ui) has a F-convergent subsequence Fg, ^{-.^Ui), 
Fr 1 (•,t^2) has a F-convergent subsequence Fg, n(-,U2), 

k k 

Ff, 2 (■) f^s) has a F-convergent subsequence F^ ^{-.^Uj,), 



Now setting n^ = n^, we see that F^^^ (•, U) F-converges to a F{-, U) for every U ^ Ai. In 
the following we will only consider the sequence e„j. and, in particular, define F' and F" 
accordingly. Furthermore, we define F{y, U) := F'{y, U) for every y and U. 
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For W CC U C il.,hy compactness of W, we always find V & Ai such that W CV CC 
U. Hence, by Lemma 13.71 we have 

F"{y, U) = sup{F"{y, V):VccU,Ve Ai} 

for every U G A{^) and y G W^''P{^,M.'^). Using, that F'{y, •) is an increasing set function, 
we can calculate 

sup{F'(y, V):VccU,Ve Ai} < F'{y, U) 

<F"{y,U) 
= sup{F"(y, V):VccU,Ve Ai}. 

But the first and the last term are equal, thus F'(y, U) = F"{y, U) = F{y, U), whenever 

The next step is to get an integral representation by showing that F, restricted to 
W^'P{Q,;M.'^), satisfies the conditions (i)-(vi) in Theorem 12.51 We immediately see the 
locality (i), by Lemma 13.81 and the growth condition (iii), by Lemma 13.51 Furthermore, 
since the F^^ are translation invariant in y, so is F, which yields (iv). To get the lower 
semicontinuity (v), just remember that weak convergence in W^'P{Q,,M.'^) implies strong 
convergence in LP(ri,R ) and that r(X)-limits are sequentially lower semicontinuous with 
respect to the convergence in X. 

To get the measure property (ii), it is enough to show that we can apply the De- 
Giorgi-Letta criterion (Lemma 12. 6p with p = F(y, •). Obviously F{y, •) is an increasing set 
function and F{y^%) = 0. Remark that for every W dd U U V{W^ U, V open), there are 
open sets U', V such that U' CC U, V CC V and VF C [/' U V , which is easily seen by 
the compactness of W . Hence the subadditivity follows from the Lemmata 13.61 and 13.71 
The inner regularity is explicitly given by Lemma 13.71 The superadditivity we can show 
directly. Take a sequence y^ G Be^ {U U V) such that y^ — )• y in L'p{Q.;W^) and 



F(y,C/UF)= limF,„^(yfc,[/UF). 

fc— >-CXD '' 



Then, 



F[y, UUV)> liminf F, (y^, U) + liminf F, (y^, V) 

k—)-oo " k—^oo " 

>F{y,U)+F{y,V), 

since U CiV = 0. Hence, we can apply the De-Giorgi-Letta criterion and obtain (ii). 
Finally, condition (vi) states that for every M G M*^^", z, z' G ^ and r > such that 
Br{z),Br{z') C ri, we have 

F{yM,Br{z))=F{yM,Br{z')), 

if we set yM{x) = Mx. By inner regularity, it is enough to show that, for any r' < r, 

F{yM,Briz))>FiyM,Br'{z')). 

Let T/fc G I3e„ {Br{z)) such that y^ —^ yM in LP(r2;M"') and 

lim F, {yk,Br{z)) = F{yM,Br{z)). 

fc— >oo '^ 

Denote by a^ the only point in C^^ (1 Qs„ {z' — z). Then define 
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,, \yk{x-ak) + Mak if x G (S^/(z')f "'^ 

Uk[X) = < 

I Mx else. 

If k is large enough, then x — aj. ^ {Br{z))e„ , Uk G B^^ {Bj.i{z')) and 

Viifc(x) = Vyk{x - ak) 
for all X G {Bri{z')Y'^k . Hence, 

Furthermore, we have Mak — )• M{z' — z) and yk{' — o.k) — ^ M(- — (z' — z)) in L^{Bf.i{z')] 
and therefore u^ — )• yM in 1/^(^2; M). Hence, we get 

F{yM,Br'{z')) < liminfFj^ {uk,Br'{z')) 

fc— >-oo '^ 

<liminfF£„ (yfc,5r(2)) 

fc— >-oo *" 

= F(yA/,Sr(^)), 

and (vi) is proven. 

Consequently, we can apply Theorem 12. 51 to the restriction of F to W^'P{Q; W^) x A{Q). 
In particular, there is a continuous function /: M — ;• [0, oo) such that 

F{y, U) = j f{yy{x)) dx if y G W^^^iSl; M!^) 
u 
and 

< /(M) < C(l + |M|P) for every M G M'^''^. (3.7) 

The asserted lower bound on / is instantly obtained, if we apply Lemma [3. 4 1 to yM and use 
the integral representation. And finally, / is quasiconvex by Theorem 12. 4[ since F(-,0) is 
sequentially lower semicontinuous with respect to weak convergence in Vl^^'^(0;M ). 

Now, let U have Lipschitz boundary. Take y G LP(ri;M'^) n Tyi'P(C/,M'^). By Lemma 
13.71 we have 

F"{y, U) = sup{F"(y, V):VC(lU,V^ Ai}. 

Using that F'[y, •) is an increasing set function, we can calculate 

sup{F'(y, F) : y CC [/, F G ^i} < F'{y, U) 

<F"iy,U) 
= sup{F"(y, V):VccU,Ve Ai}. 

But the first and the last term are equal, thus F'{y,U) = F"(y,U) = F(y,U). If y G 
LP{n;m'^)\W^'P{U,R'^), then oo = F'{y,U) = F"{y,U) = F{y,U) by Lemma [331 Hence, 
Fen (■'^) r(L^(r2;M ))-converges to F{-,U). To get the integral representation for y G 
LP^Q-jW^) n W^'P{U,M.'^), observe, that since U has Lipschitz boundary, we can find a 
function v G W^^p{Q; R'^) such that 

y{x) = v{x) for almost every x € U. 

Then, by Lemma 13.8^ 



F{y, U) = F{v, U) = J f{Vv{x)) dx = J f{Vy{x)) dx. 
u u 

D 
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4 The boundary value problem 

While loading terms can be included in our results so far without difficulties, the restriction 
to deformations with preassigned boundary values is more subtle. 

4.1 Statement of representation result with boundary conditions 

Suppose g S 1^^'°°(]R'^;IR'^) is a boundary datum. We will then always choose the precise 
representative for g and thus assume that g is continuous. We define the admissible lattice 
deformations Bs(U,g) as the functions in Bs{U), that satisfy the boundary condition 

y{x) = g{x), whenever x € dCeiU). 

The correspondingly restricted discrete functional is 

I oo otherwise. 

Assume that £„,, and / are as in Theorem 13. 1^ let us for simplicity write just £k in the 
following and set 

I oo otherwise. 

In analogy to Theorem 13.11 we then have: 



are 



Theorem 4.1. Suppose Assumptions\^ and\^ are true, g G VF^'°°(R'^;M'^) and F^ , F|J, 
as above. Then Fi^{-, U) T{LP{n;R'^))-converges to F9{-, U) for every U G Al{^). 

4.2 Improved estimates on interpolations 

We start by improving Proposition 13.31 for sequences in Be{U,g). This is possible, because 
now we can control what happens near the boundary. Note, that now we can naturally 
define the interpolation y on all of [/, namely, we just extend y by the discretization of g 
before we interpolate. 

Proposition 4.2. Let U G Aiift) and y^ G Be^{U,g). Then y^ ^ y in LP{U;R'^) if and 
only ifyk^y in L^(C/;M'^). 

Proof. First, let y^ — > y in LP(C/;M ). Choose some open bounded set U' with Lipschitz 
boundary and U CC U' . Extend the functions by defining yk{x) := g{x) and y{x) := g{x) 
for X G U'\U . So yk G Bs^{U' ,g) and, since g is Lipschitz, we have yk ^ y in L^([/';M ). 
But then by Proposition 13.31 we get y^ — t- y in LP{U;W^). 

Now, let jjk ^ y in LP{U;M.'^). Let Aj : M*^ — )• [0,1] again denote the cell-periodic 
functions such that, with zi = A{—^, . . . , — i), 

i{Zi - zi)). 



yn{x) =^^i[ — ] yn{x + £«( 



Of course, Aj > and the Aj add up to 1 in any point. Define 

Wn,i = \x eU: \i( — \ > - and Aj f — j < a iov j y^ i 
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where a will be chosen suitably later, and note that, for every x € Us„ , the ratio - — rg — m] 
is independent of n and x and positive since x G <5e„, x —^ x -\- Zi implies that Aj( — ) — ?> 1 
and Aj( — ) — ?• for j ^ i- Next, extend y and yk by g as above and define 



Pny{x):= j y{h) 



_ ,dh. 

Of course, we have \\Pny — y\\LP(u-M.'i) ~^ ^- Hence, it suffices to show \\Pny — UnWLpm-^ 
0. For X € Wn^i we have 

|y„(x - en{Zi - Zi)) - Pny{x)\ >- \yn{x) - Pny{x)\ 



X] ^J \ ) 1^"*^^ ~ ^"*^^* ~ ^j)) ~ Pnyix) 



Since y„ and P„y ^-re constant on every cell, we thus have 

1|| P II - 1 I^^J^ II p II 

l^-^En ' ' '' n,l\ 

\UeJp 
— ~ 1 \\yn{- — £n\Zi - Zi)) — Pny\\LP(Ue„nWni;M.'i) 

\Ue„nWn,i\p 

+ '^2^ Wy^iV ~ ^n{Zi — Zj)) — ^«?/|lx,p((7^^;]Kd) • 
1 

But — —^ — T > is independent of n and 

\u,„nw„,^\p 

\\yn{- — £n{Zi — Zi)) — Pny\\LP(Ue„nW„x,^'i) 

< \\y — Pny\\LP{U;R'') + WVni- — £n{zi — Zi)) — ?/||/,p({/;Kd) 

converges to 0. To control the remaining sum, we estimate 

\\ynV ~ £n{Zi — Zj)) — Pny\\ip(^if^^-^d'^ 

< \\yn{- — £n{Zi — Zj)) — Pny{- — £n{Zi — ^i)) IIlp((/^^ ;Rd) 
+ \\Pny[' ~ £n[Zi — Zj)) — Pny\\i^p(^u^^.^^d^ , 

where the second term goes to and the first term is estimated by 

\\ynV — £n{Zi — Zj)) — Pny{- — £n{Zi — ^jj) IIlp{(7e„;R<*) — 11^™ ~ -fri?/llLP((7^n;R'*) • 

Altogether we obtain 

\\yn — Pny\\LP{U;R'^) 

< \\yn - Pny\\LP{U^n\U^^;R<i) + \\yn " -Pn?/ |lLP((7^^;]Rd) 

Ilp(C/="\C/£„;R'*) + ^^(^ ~ ^) 11?^" ~ -" ny\\LP{Ue„] 



= (1 + 20(2"* - 1)) ||y„ - Pny\\LP(U^n\u, -Rd) + 20(2^^ - 1) ||y„ - Pny\\LP(U, ■R'*) + ^(1) 
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As near the boundary we can calculate 

1 
WVn - -Pny|lLP([/£n\(7^^;Rd) < |C/'^"\f/e,J P \\g\\^ + \\Pny\\LP{U^'^\U,„; 

for a = 2dqn" we finally get 



D 



Remark 4.3. The proof shows that without boundary conditions, i.e., for a general sequence 
Vk € Bs^^iU), U e A{^), we stih have y^ ^ y in Lf^^{U;R'^) if and only if y^ ^ ?/ in 

4.3 Proof of the boundary value representation result 

Proof of Theorem \4-l\ Fix U € ^i(ri). We start with the lim inf-inequality. Let yk,y ^ 
L^{Q,; M ) such that yk ^- V- We can assume that 

liminfF/ (yfc,f/) < 00, 

because otherwise there is nothing to show. For some subsequence we then get 

hminf F4(yfc,f/) = lim F| (y,,,C/). 

But since F^^, < Fi^ , we can argue as in Lemma 13.41 to see that y € W^'^{U]Mr') and, 

for any V <Z<Z U, that yk^ -^ y in W^'^{V;M.'^). Using Proposition 14.21 we see that ykj 
converges strongly in LP{U;W^) and, since Vy^^ is now bounded in LP{U;W^), weakly in 
W^'^{U; M ) to y. Regarding the boundary condition, there are open neighborhoods Vi of 
dU, where jjki is an affine interpolation of g. Namely, V/ is the interior of the union of all 
cells g^^^ , with Q^ n aC/ / 0. Then 

sup lijkiix) - g{x)\ < sup \yki{x) - g{x)\ < Ceki 
xedu xeVi 

since g is Lipschitz. Denoting the trace operator by T, we thus have Ty^^ ^ Ty = Tg in 
LP(aC/;M^) and hence y € g + Wq'^{U;M'^). But then, we can calculate 

F^(y,C/) = F(y,[/) < hminf Fe,(yfc,f/) < hminf F|^(yfc,C/), 

and have indeed proven the lim inf-inequality. 

To get the F-convergence result, we now proof the lim sup-inequality. Let us first assume 
y{x) = g{x) + Tpix), for every x G U and some Tp G C^(C/;M ). Then 

F3{y,U)=F{y,U)<cx^. 

So, there exists a sequence Uk S Bsf.{U) such that n^ — )> y in LP{Q.]W^) and 

limF,„(nfc,C/) = F(y,[/). 

k^oo 
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Let 6 > 0, and then choose U' such that supp'i/' <Z U' CC U and |C/\[/'| < 6. We now use 
a cut-off argument similarly as in the proof of Lemma 13.61 Fix A^ € N and define 



U, = [.eU:aisti.,U')<^:^^^^] 



Then choose the cut-off functions (pj E C^([/j+i; [0, 1]) with ipj = 1 on Uj and ||Vv3j|| < 
CN and set 

gkix) = g{a), if a G Qski^) n Ce^ and 

( ^_\^A^)^k{x) + {l-^j{x))gk{x), ifxGC/^*, 
|iifc(j;) otherwise. 

As in the proof of Lemma [3.61 we calculate 

Fe,{wk,j,U)<Fe,{uk,U) + C{\\Vg\\l^ + l)\U\U'\+ei Y. ^ceii(Vt/;fc,,(s)), 

'x(i(c'xwj)r 



^j,k 



with Wj = Uj-\-2\Uj-i, estimate Sj^k by averaging, choose j{k) suitably and thus get 
limsupF,,(wfc^^(fe),C/) < F{y,U) + C5+j^—-^. 

Since we choose j[h) < N — 3, we have Wkj(k) ^ ^Ski^^d) ^o^' ^^^Y ^ large enough. Fur- 
thermore, Wkjik) ~^ y hi L'^{Q;M. ) since ^p has support in U' . Hence, 

r-limsupF/Jy, U) < F'^{y, U) + 6C + -^. 

fc— >CJO -'* "^ 

Let (^ ^ and N ^ oo. 

In the general case y\u G 5 -|- VFq'^([/;]R'^), take y; such that yi\u G g -|- C^([/;M'^) and 
yz ^ y in W^^p{U;R'^) and in LP([7;M'^). We get 

F- lim sup F/^ iy,U) < lim inf (F- lim sup F/^ (y; , U) ) 

fc— >oo (— ^00 fc— >oo 

<liminfF3(yz,C/) 

i— )-oo 

= F5(y,[/) 

by the lower semicontinuity of the F-limsup with respect to LP(f];]R )-convergence and 
continuity of F^^-, U) with respect to VF^'P(f7;R'^)- convergence. D 

4.4 The limiting minimum problem 

The following theorem is important in two ways. On the one hand we gain insight into the 
F-convergence result, on the other hand we will directly need it to get the homogenization 
result in Section [5l 

Theorem 4.4. Under the assumptions of Theoreni \4-l\ we have 

mmF3{y,U)= lim (inf F/ fy, f/)). 
y k^foo y 
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Furthermore, any sequence yj. with equibounded energy is pre-compact in LP{U;M. ) and if 
we have a sequence satisfying 

lim (infF|Jy,C/)) = lim F|^, (?/fc, C/) , 

then every limit of a converging subsequence is a minimizer of F^{-, U). 

Proof. Fix g(, U and write Gk{y) = Fikiv^ U)^ G{y) = F^{y, U). Let yt be a sequence with 
equibounded energy Gk{yk)- By Assumption [1] and Proposition 13.21 we obtain that 

f IVm]" dx < c. 

Furthermore, using the boundary condition, we have 



J\Vykfdx<C. 



u 
A Poincare-type inequahty involving the trace yields 

<C + C\\g\\^n'^~^{dU)^ <C 

and so y^i —^yin LP[U]W^) for some subsequence ki and some y G W^''p{U]W^). Then, 
by Proposition [Ol yk, ^ y in LP{U;R'^). 

Now from Theorem 14.11 we infer that Gk r(LP(Q;IR'^))-converges to G. But then G^ 
also r(L*'(J7;IR'^))-converges to G. Here the existence of recovery sequences is immediate 
as L'P{U;'R. )- implies L^(C/;M )-convergence. As for the lim inf-inequality, if y^ ^- y in 
LP{U]M.'^), where the energies Gk{yk) are without loss of generality assumed to be equi- 
bounded and, in particular, in y^ S Bef.{U,g), we can extend the functions by defining 
yk{x) := g{x) and y{x) := g{x) for x S f]\[/ without changing their respective energies. 
Since then y^ ^ y in L*'(il;M ), we have indeed that liminik^oo Gk{yk) ^ G{y). The 
remaining part of the proof now directly follows from Theorem 12.31 D 

5 Proof of the main results 

5.1 The F-convergence results 

To simplify notations, we define Ph{x) = x-|- A(0, h)'^ and P/, = Ph,(0). First, we will prove 
the following lemma. 

Lemma 5.1. The limit 



exists for every M € W^^'^ 



lirn^ ^ inf <j Y^ VFccii(Vy(x)) : y G Bi{Pn, yhi] 
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Proof. Let us define 



G{y,U)= J2 ^ceii(V2/(x)) and 

xGic[iu)r 

fk{M) = -1 inf {G{y, Pk):ye Bi{Pk,yM)} . 
Fix M £ M''^'^ and let k,n G'N with k > n. Clioose u„ G Bi{Pn,yM) sucli that 

^G{Vn,Pn)<fn{M) + l 



n 



n 



Now, we can define 

Vn{x — na) + nMa if x G Pn{na) for some a ^ A iO,l, . . . , 



Uk{x) 



Mx 



otherwise. 



Since Vn satisfies the boundary condition, u^ is constant on every ceh. Moreover, u^ € 
l3i{Pk,yAi) and we can estimate 



fk{M)<^G{uk,Pk) 



< 



k'^ 



G{Vn,Pn) + c{\M\P+l){#{£[{Pk)r 



#{C[{Pn{an))y 



^>--.)+^^^s^ 



P.rAi 



n" 



Pi 



+ 



n'^ - (n - 2)' 



^ /«(-)- ^ ^^4^ '^-'" 



..Mi-(i-^) 



</„(M) + i+c(|Mr + l) l-(l-f)Vl-(l-l 



Thus, for every n € N, 



hmsup /fc(M) < /„(M) + - + c(|Mf + 1) ( 1 - ("l - -^ 



hence, 



Hmsup /fc(M) < hminf /„(M). 



Now, we can prove our first main theorem. 
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Proof of Theorem[Tj\ We will first show that Fe(-,0) r(LP(ri;M'^))-converges to F. Ac- 
cording to Lemma 15.11 Wcont is well-defined. By the Urysohn property of F-convergence 
in Proposition 12. 2^ it is enough to show that, for any sequence e^ — >■ 0, the function / of 
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Theorem 13.11 equals VFcont- Fix such a sequence, the subsequence e^ and the associated /. 
Since / is quasiconvex, we have for every M € M and U G Al{^) 



/(M) = y^min< 



1 



mm 



u 



If we restrict yM to a ball that contains some neighborhood of Q, we can extend it to 
a function in W^'°°{M.'^;M.'^) fl C(M'^;M'^), so yM is admissible as a boundary condition in 
Theorem 14.11 and we get the T-convergence result with boundary condition. Hence by 
Theorem I331 for /iq > and xq G M'^ such that Pho{xo) CC $7 



f{M) 



1 



1 



lim {mi{Fs^^{y,Pho{xo))- y G SejP/j,j(xo),yA/)}) 



\Pho{xo)\ fc^oo 



lim inf < e^. 



E 



W^ceii(Vy(x)): y G Be^{Pho{xQ),yM) 



xG(£^JPho(a;o)))° 

It is easy to see, that we can always find /i^ > and x^ G £e^, such that 



^/ife(a:^fe) 



U Q^.(^ 



We then know Ph^{xk) C for all k large enough, |xo — Xk\ < diamQgj, = e^diamQi, 
h{) < hk < Hq + 2ek and, that there are Nj^. G N satisfying /i^ = Nj^Ej.. Furthermore, 

C',^{Ph,{x^)) = C',^{Pf,^{xk)) and 
(4^(P,„(xo)))° = (4,(^h..(^fc)))°- 

Hence, Bs^{Phoixo) , yhi) and Bs^{Ph^ixk),yAi) are equal up to extending the functions in 
B£^{PhQ{xQ) , yM) constant on cells that intersect Ph^{xk) \ P/iq(xo). It follows that 



f{M) 



1 ,- ^ K- , 

hm — r — T mf 



I Pi I k^'^oN^h^ 



7-— 7 lim — J inf 

I Pi I k-^00 Ni 



E 



VFccii(V2/(x)): y G Be^{Ph,(.Xk),yM) 



_xe(£^JP^Jxfe)))° 



J] H^eeii(Vy(x)):yGS,,(P,,,yM)L 



,a;G{£',(Ph,))° 



where we used, that y G Be^{Phk{xk)-,yM), if and only if y(-+a;fc)-Mxfc G Bef^iPhk,yM) and 
that the discrete gradient of y at a point x equals the discrete gradient of y(- + Xfc) — Mxk 
at x-Xk- In asimilar way y G Be,^{Phk,yM) if and only if y' G Bi{PN^,yM) and Vy'{x) = 
Vy{ekx), where y'{x) = -^y{ekx). Hence, 



/(M) 



Efe' 



lim — T inf 



I Pi I k'^'ioNi 



Wcont(M). 



Y, W^ceii(Vy(x)) : y G ei(P^, , ^m) 
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In order to prove that also F^{-,Q,) r(L^^^(0;]R'^)/R)-converges to F, we only need 
to verify the lim inf-inequality as the existence of recovery sequences immediately follows 
from the first part of the proof since convergence in LP[Q;W^) implies convergence in 
Lf^^{n;W^)/R. But if e„ ^ and y^^ ^ y in Lf^^(ri;M'^)/R, then there exist c^ G M 
such that, for every U € Al{^) with U CC fi, y^^^ — Cn ^ y in LP(C/;M'^), so that by the 
previous result 

limmiFe„{ye„,n) = lim inf F^^ (?/e„ - Cn,^) > lim inf F^^ (?/£„ - c„,C/) > F{y,U). 

n— >oo n— >-cxD n— >-oo 

Without loss of generality we may assume that lim inf F^^, (y^^, , fi) < oo. Since for any 

V e A{il.) with y CC il there exists U € Al{^) with V C U CC il., we then deduce from 
Lemma 13.41 that y E W^'''^{V;W^) with ||2/|[jyi,pcy.|gd) bounded uniformly in V £ A with 

V CC il, hence y C VK^'^(ri;R ). Then invoking Lemma [3.71 and passing to the supremum 
over U G ALi^) in the above inequality yields 

limmfF,^{y,^,n)>F{y,n). 

fe— ^oo 

D 

Proof of Theorem \1.3[ Theorem 11.31 is a direct consequence of Theorem 14.11 and Theorem 
11.11 where the limiting energy density / has been identified as Wcont- D 

Proof of Theorem ] 1. Si Suppose y^ is a sequence with equibounded energies F^^{yk). By 
Proposition 13.21 and the growth assumptions on Wceii, for every U G A{^) with U CC O 
we have 

/ \Vykf<CF,^{y,) + C\n\ 
Ju 

uniformly bounded for sufficiently large k. Choose Uq € ALi^) connected and with ^ 
Uq CC ^. As C/q is connected, by Poincare's inequality we find Cfc € M such that y^ — Ck is 
pre-compact in LP{Uo;M.'^). But then indeed for any connected U G Al with Uq C U CC ^ 
the Poincare inequality 

WVk - Cfc|lM/i,p(c/;Kd) < C ||Vyfc||^p(f/.igd) + \\yk - Ck\\LP(Uo;M'i) 

yields that y^ — Ck is pre-compact in L^(C/;R ). Exhausting Q, with a countable number 
of such domains and passing to a diagonal sequence, we find a subsequence y/j,^ such that 
yk„—Ck„ converges in L^^^(il; M ). By Remark 14. 3 1 we finally obtain that y^^ —Ck„ converges 
inLf^^(0;M'^). D 

Proof of Theorem \1.4\ This is immediate from Theorem 14.41 D 



Proof of Corollary IJ.51 This is a direct consequence of Theorems II. Ij II. 2| II. 3| 11.41 and 
[2:31 D 



Proof of Theorem \1.8[ If in addition to Assumptions [T] and [2] Assumption |3] holds true, 
we can apply |CDKM06| Theorem 4.2] with A = {C[{PnJ)°. It is easy to see that the 
boundary of A as defined in |CDKM06] equals 5£i(P/v^) U £i\£'^(PjVj.), but of course the 
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second part does not change anything. This shows that there is a neighborhood U of 
SO{d), such that for every M €U 






Y. H^ceii(Vy(x)): y € ^i(Pjv,,yM) 

^ WcelliMZ) 



|det^, 
Wcb{M). 



D 



5.2 Approximation of general continuum densities 

Next we prove Propositions 11.91 and ll.lOl 



Proof of Proposition IJ.91 At variance with our previous decomposition procedure, we now 
choose any simphcial decomposition S of the cell yl[— ^, ^)'^ into d-simplices all of whose 
corners lie in A{— 2, 2} • For F = (/i, . . . , /2d) € M we then interpolate the mapping 



^<!-^,U ^IR^ xi^fi 



1 i''^ 

2'2 



afiine on each simplex in order to obtain 

Up : A 
Then Wceii is defined by 

VFeeii(F):= / ViVuF)dx. 

As every corner Zj^ , . . . , Zj^ of 5 € 5 lies in A{—^, i}'^, we have 



on S. Thus, lli^^ll = max \SJuf{x)\ is a norm on Vq and we calculate 

x&A[-UY 



^cell(i^)> / c\VuF\''-ddx 



>c||Ff -c' 
>c\Ff-d, 

and on the other hand 



29 



This means Wcdi satisfies Assumptions [T] and [2j From Theorem 1 1.1 1 we then deduce that 
^^cont(M) = -^ Jirn^ W'""^] ^ ^cell(Vy(x)) : y € ^i(Piv, ^A/) 



= -4- lim -^ \detiA)\ {N - 2)''V{M) 

|det j4| n^oo Jy"- 

= V{M) 

due to the quasiconvexity of V. D 

Proof of Proposition \1.1(A Any F € Vq can be decomposed orthogonaUy as i*" = F'Z + F" 
with unique F' G M^^"' and F" G (R^'^'^Z)-^. Set 

Ty,eii(F) = |det A\ Q {J{F'YF' - ld\ + \F"\^ + x(F), 

where % is any frame indifferent function satisfying Assumptions [1] and [2] with p > d which 
is non- negative, vanishes near SO{d) and is bounded from below by a positive constant on 
0{d) \SO{d), 0{d) = 0{d)Z. Then also VFceii satisfies Assumptions [1] and [2] with the same 
p. Noting that, for M £ W^""^, (MF)' = MF' and {MF)" = MF", it is not hard to verify 
that Wceii also satisfies Assumption [3] with 



D'W,,ii{Z){F, F) = 2 |det A\ Q ^^— + 2|F 



But then 



1^9... ..,w,. ,.. 1 ^2... ,^w,.^ ,.^. ^/Af^+Af 



-D'WcB{ld)iM,M) = -^^-^^D'W,,ii{Z)iMZ,MZ) = Q ( \ = Q{M). 

D 

5.3 Extension to finite-range energies 

We briefly comment on more general long-range interactions. Suppose A = {zi, . . . , ^2^, . . . , zat} C 
C is any fixed finite set, where zi, . . . , z-^d still denote A^ — i^, 2}'^- For y € ^£(^) we define 
yi = y{x + ezi). With x and y as before, i.e., only depending on yi, . . . j/2d, let now 

The lattice interior {C'^{U))° and boundary dC'^{U) now have to be shrunk respectively 
enlarged to a whole boundary layer, according to the maximal interaction length in A. 
Assumptions [1] and [2] are then replaced by the estimate 

C iFf - c' < t^supcr-ccll(i^) < C"(JF|^ + 1) 

for constants c,c',c" > and aU F € M''^^ which satisfy F' e Vq, where F' e M'^^^'' 
denotes the left d x 2 submatrix of F. Note that the lower bound in particular allows 
for arbitrarily weak long range interactions. As the interpolation we used only depends on 
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the d X 2'^ values of the corresponding lattice cell, this implies that we get the standard 
estimates for the gradients in Proposition 13.21 only on this part of the discrete gradient. 

It is important that the interaction range is bounded by Ce, so that, e.g.. Lemma 
13.51 and its proof still work. In the estimates of the error Sj^n in, e.g.. Lemma 13.61 it is 
important that £~^ \u{x + ezj) — u{x)\ and thus the discrete gradient can be bounded by 
a fixed finite sum of smaller d x 2 discrete gradients of some cells near x. Hence, we still 
have the estimate 

si Yj |Vn(x)|^ <C f \Vu{x)\P dx 

Note that according to our enlarging of the lattice boundaries, also the cell formula for 
the limit density will now involve a sequence of minimizing problems with affine boundary 
conditions on a boundary layer. 

We finally remark that the statement on the applicability of the Cauchy-Born rule 
translates naturally, as the main ingredient does, see |CDKM06l Theorem 5.1]. 

5.4 Extension to multi-lattices 

It is also possible to generalize these results to certain non-Bravais lattices, namely to 
multi-lattices of the form £ U (si + £) U • • • U {sm + ^)) in the following way: We still 
consider C to be our main lattice. But now we have m additional atoms in each cell, which 
we describe by the 'internal variable' s{x) € R , such that es.j describes the distance of 
the j'-th atom to the midpoint of the cell. Of course s can be identified with a function, 
that is constant on every interior cell and is outside and thus lies in some L'?(ri;M ), 
1 < (/ < oo. The new cell energy depends on md additional variables and we now consider 
the growth condition 

C(|M'|^ + |sn -d < I^supcr-cell(M, s) < c"{\M\P + |s|'' + 1) 

for M € R and s € R . It is now natural to have a F-convergence result with respect 
to strong- LP-convergence in the first and weak-L'?-convergence in the second component. 
As we will see in a moment, it turns out that we have to consider a combined boundary 
value and mean value problem. For this we define Be{U, g, sq) to consist of all pairs (y, s), 
such that y G Bs{U,g) and s € L''(0;R'^^'") is constant on every interior cell of [/, is 
outside and has mean value sq on the union of interior cells of U . 

In analogy to Theorem 11.11 we now have Theorem 11.61 The proof of this theorem is 
similar to the proof of Theorem ll.il But there are several things that need to be addressed: 

First of all, the weak topology on L'^ is not given by a metric. But, as discussed in 
|DM93| in detail, this is not a big problem, since our functionals are equicoercive and the 
dual of L^ is separable. In particular, we can describe F-convergence by sequences and the 
compactness and the Urysohn property are still true. Next, we need an advanced version 
of our integral representation result: 

Theorem 5.2. Letl<p,q<oo and let F : W^'P{n;R'^) x L'?(0;R'^^"") x ^(0) -^ [0,oo] 
satisfy the following conditions: 

(i) (locality) F{y,s,U) = F{v,t,U), if y{x) = v{x) and s{x) = t{x) for a.e. x € U ; 
(ii) (measure property) F{y, s, ■) is the restriction of a Borel measure to A{Q); 
(Hi) (growth condition) there exists c > such that 

F{y, s,U)<c f \Vy{x)f + |s|'? + 1 dx; 



u 
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(iv) (translation invariance in y) F{y, s, U) = F[y + a, s, U) for every a G 



pd . 



(v) (lower semicontinuity) F{-,-,U) is sequentially lower semicontinuous with respect to 
weak convergence in W^'^{Q.]M. ) in the first and weak convergence in L'^(rj;M ) 
in the second component; 

(vi) (translation invariance in x) With yuix) = Mx and s{x) = sq we have 

F{yM, s, Br{x)) = F{yM, s, Br{x')) 
for every M G W^'^'^, sq G M'^''"', x,x' (^Q. and r > Q such that Br{x),Br{x') C 0. 
Then there exists a continuous f : M'^^'^ x M'^^™ — > [0, oo) such that 

0<f{M,s) <C(l + |M|P + |sH 
for every M G R'^'^'^, s G W^"^"^ and 



F{y,s,U)= I f{Vy{x),s{x))dx 
u 

for every y G W^'P{n;R'^), s G L9(J];M'^><™) and U G A{n). 

The proof in |BD98| for the pure Sobolev version of this theorem readily apphes to this 
more general statement. (Note that continuity of / then follows from seperate convexity.) 
Most of the lemmata then translate naturally. We just want to comment on some details 
in Lemma 13.61 The recovery sequences now contain additionally some tn ^^ s, rn ^- s 
corresponding to U,V respectively. We define 

qn,j = XUj{x)tn(x) + (1 - XUjix))rn{x), 

and then choose j{n) as before to define s„ = qn,j{n)- The only part that is not immediately 
clear now, is the convergence s„ ^ s. To prove this, let 93 G L"? (fi;R'^^'"). We now split 
(p into several parts we can control 

N-l 

c/P = Vn + ^XU' + ^Xn\u^ + Yl '^X(c/,+i\I7-),„ • 

i=o 

Here the ipn contain all the remaining parts. We see that Vn — ^ strongly in L"? as long 
as \dUj\ = for every j, so this is true up to changing the sets Uj a little bit. But then we 
also have 

strongly in L'^ . The advantage is now that on each set {Uj+i\Uj)s„ we have either Sn = tn 
or Sn = Tn, possibly changing with n. But in both cases we have weak convergence to s, 
hence 

Sn{x)ip{x)x^jj^^^\Tj-)^^{x)dx ^ I s{x)ip{x)xu,+^\u,{x) dx. 
n n 

And, putting it all together, we get 

Sn{x)ip{x) dx ^ / s{x)(p{x) dx. 
n n 
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Another important step is to adjust Theorem 14. II and Theorem 14. 4| so that additionally 
to the boundary values for y we have a fixed mean value for s, i.e., we consider B^{U, g, sq) 
instead of BgiUjg) in the discrete setting and add the constraint 



s{x) dx = So 
u 

in the continuum setting. To get the lim inf-inequality just notice that for Sk -^ s with 
{yk,Sk) G Be{U,g,so) we have 



/s{x) dx = lim -r Sfc(x) dx 
u 

lim — ^ f Sk{x)dx 

k—^oo \U \ J 



U U 



Ue, 



= hm -ttt^sq = Sq. 
k—^oo \U I 

The lim sup-inequality is a little more subtle. We have a function s G L'?(ri;R ) with 

s{x) dx = Sq 



u 
and a recovery sequence without this mean value Sk —^ s. Let us write 



-h Sk{x)dx + ik = So, 



so that ^fc — 7> 0. We now adjust the Sk adequately. Define 

tk{x) = Sk{x) + ^fc-r^XVfc- 

If Vk is a union of cells with some distance to the boundary of U , then, for k large enough, 
the tk are admissible functions and do not interact with the adjustments on y. We have to 
make sure, that tk -^ s and 

Imisn-^ Fe^^{uk,tk,U) <\vms\\Y) Fe^{uk, Sk,U). 

k—^oo fc— >oo 

The weak convergence is true, if |Vfc| — ?• and |Vfc| > c^k for some c > 0. For the second 
estimate, we have to choose the Vk a little more carefully to avoid concentration of the 
energy. Choose sequences 77^ — > 0, Lfc — )• oo such that rjk > c^fc, -^ — ?> co and LkTjk — ?• 0. 

Then take L^ G N disjoints sets Wk^i C U, that are unions of cells, such that |M/fc,z| is 
independent of / and is roughly equal to %, which means 

crjk < \Wk,i\ < Cr]k, 
with C, c > independent of k and /. This is possible as % — > 00 and Lfc??fc — ^ 0. Then, 



fc 
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we can choose l{k) and set V^ = ^fc,i(fe)) such that 

Wco\\{^Uk{x), Sk{x)) + Wccii{Vuk{x),Sk{x) + ^Trfy) dx 

I ^'' f - - \U \ 

<-^^ / Wcell{S/Uk{x),Sk{x))+Wcen{S/Uk{x),Sk{x)+^kT^^)dx 



'=V,, 



due to the growth condition. So the error goes to zero with L^ —^ oo. The rest of the proof 
translates naturahy. The most important observation is the equahty 

f{M,so) = i^minj j f{Vy{x),s{x))dx: y - yu G <'^(C/;M''), 

U ^ 

= -^ min I F{y, s,U):y- yM G <'^(C/; K'^), 

s G L'?(C/; M'^^™), / s{x) dx = soV 

which is of course a consequence of the lower semicontinuity properties. 

Proof of Theorem \1.7[ Fix y € VF^'P(0;IR'^) and without loss of generality fix a version of 
y that is finite everywhere. 

Due to the growth condition and the continuity, we know that the infimum in 

inf Wcont{M,s) 

is actually a minimum for arbitrary M and that the function 

M^ min Wcont (M,s) 

is continuous. Obviously, we always have the inequality 

WcontC^y{x),s{x))dx > min WcontC^y{x),s)dx. 

J J s^Wx^ 

n n 

We now want to show, that there always exists an L'^-function s where this is an equality. 
The idea is of course to choose s{x) as a minimizer of s i-^ ^cont(Vy(x), s). The key point 
is to ensure measurability. We will do this by using the theory of measurable multifunctions 
as developed, e.g., in |FL07| . Define 

e(M) = {se M'^^'": WcontiM, s) = min W,ont{M,t)} 



and set T{x) = 0(Vy(x)). Due to the continuity and the growth of Wcont, the set T{x) 
is closed and non-empty for every x G fi, hence F: fi ^ 7-'(M ) is a closed- valued 
multifunction. 
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Next, we want to show that T is measurable, in the sense that 

r-(C7) = {xen: r(x)nc/0} 

is Lebesgue-measurable for every closed set C C M ™. To this end, we will first show 
that 0"(C) is closed. Let M„ G ©"(C), M„ -^ M and choose s„ G ©(M„) n C. Using 
the growth of Wcont and since the M„ are bounded, the Sn are also bounded. So for some 
subsequence we have Sn^. -^ s and s a C. Furthermore, 

WcontiM,s)= lim Wcont{Mn„SnJ 

fc— )-oo 

= lim min Wcont{Mni^,t) 



min Wcont{M,t). 



This proves s G G(M) nC, so ©"(C) is closed and r-(C) = (Vy)-i(G-(C)) is Lebesgue- 
measurable. Now, we can apply |FL07| Thm. 6.5], to get a measurable s: ^ R'^^™, 
with 

VFcont(Vy(2;),s(x)) = min Wcont(Vy(x),t) 



and s G L'?(J];R'^^™), since 

/ |s(x)|'' dx<C I W^cont(Vy(x), s(x)) + 1 (ix 
n n 

= C min Wcont(Vy(2;),s) + Idx 
n 

<C f min |Vy(x)f + lsl^ + l(ix 

J sgR'ix™ 

<C f \Vy{x)\P + ldx. 



d\ 



It remains to justify the F-convergence result for F^ ™™. Suppose yk ^- y ^ W^'^{Q,; 
strongly in LP^Q; R'^). Choose Sk G L9(Q; R'^x'") with F^^-"^"'^(yfc, f]) < ^^..(yfc, Sfe, Jl)+A;-i. 
Without loss of generality assuming that F^~^^^{yk,^) is bounded, by passing to a subse- 
quence (not relabeled) we may assume that s^ -^ sq in L^. But then Theorem 11.61 shows 
that 

lhnmfF,\:'^'^{yk,n) =\immf F,^{yk,Sk,n) > F{y,so,n) > F'~'^'^{y,n) 

fc— i-oo fc— s-oo 

by the first part of the proof. On the other hand, if y G W^''P{Q;W^) is given, choose s G 
L''(0;R'^x"') according to the first part of the proof such that F''~"^''^(y, 0) = F{y,s,n). 
Then if (y^, s^) is a recovery sequence for (y, s) from Theorem 11.61 we obtain 

limsupFi-^'^iyk,n) < limsupF,,(yfc, s^, 0) = F(y, s, fi) = F'~^'^iy,n). 

k—^oo fc— >oo 

D 
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